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Abstract 

The results of an analytical study of the nonlinear 
and buckling response characteristics of curved panels 
subjected to combined loads are presented. Aluminum 
and laminated composite panels are considered in the 
study and a flat and shallow curved panel configurations 
are considered as well. The panels are subjected to com- 
bined axial compression and transverse tension or com- 
pression loads or combined axial compression and in- 
plane shear loads. Results illustrating the effects of var- 
ious combined load states on the buckling response of 
the panels are presented. In addition, results illustrating 
the effects of laminate orthotropy and anisotropy and 
panel curvature on the panel response are presented. The 
results indicate that panel curvature can have a signifi- 
cant effect on the nonlinear and buckling behavior of the 
panels subjected to combined loads. Results are includ- 
ed that show that geometrically perfect panels do not ex- 
hibit bifurcation points for some combined loads. 
Results are also presented that show the effects of lami- 
nate orthotropy and anisotropy on the interaction of com- 
bined loads. 

Intr oduction 

To achieve high performance in the next genera- 
tion of aerospace structures, it will be necessary to de- 
velop reliable, high-fidelity design criteria for the 
buckling behavior of shell structures subjected to com- 
bined mechanical and thermal loads. For the past thirty 
years, buckling loads for shell structures have been ob- 
tained for design purposes, to a large extent, by con- 
ducting linear bifurcation analyses on idealized, 
geometrically perfect shell structures and by using em- 
pirical "knockdown" factors that account for the effects 

of initial geometric imperfections. 1 In many cases, bi- 
furcation-buckling interaction curves for shells subject- 
ed to combined loads have been approximated by a 
linear relationship between the bifurcation-buckling 
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loads for each component of the combined loads. This 
type of approximation neglects the nonlinear interac- 
tion of the combined loads and results in an overly con- 
servative prediction of the bifurcation-buckling load. 

Often, empirical knockdown factors for shells 
with complex geometry and subjected to combined 
loads are unavailable, and analysis-based predictions of 
the sensitivity of the panel response to initial geometric 
imperfections and other complicating effects must be 
used. Ideally, highly accurate predictions of the re- 
sponse can be obtained when the initial geometric im- 
perfections, materia] properties, thickness distribution, 
boundary conditions, and other features are known to a 

high level of accuracy. 2 When these quantities are not 
well known, analysis-based lower-bound estimates that 
represent the effects of initial geometric imperfections 
in nonlinear analyses must be used. A common exam- 
ple of this approach is to use a bifurcation-mode imper- 
fection shape to represent the initial geometric 
imperfections. This lower-bound, nonlinear-analysis- 
based approach has been applied to the Space Shuttle 
superlightweight external liquid-oxygen tank. 3 ' 7 This 
approach is inherently ultra-conservative and leads to 
an overweight structural design, but is often the only 
meaningful approach available to designers and ana- 
lysts. 

To complicate matters, recently it has been shown 
that geometrically perfect, shallow cylindrical panels 
that are subjected to combined axial compression along 
the generator and transverse tension or compression do 
not exhibit bifurcation points. 8 This result was ob- 
tained by using an approximate special-purpose, single- 
mode nonlinear analysis and indicates a major short- 
coming of the traditional shell stability design approach 
that is based on linear-bifurcation buckling analyses or 
results. The main objective of the present paper is to 
examine the nonlinear behavior of curved panels sub- 
jected to combined loads with a high-fidelity, nonlin- 
ear-analysis tool. This objective is achieved in the 
present paper by first describing the high-fidelity anal- 
ysis tool, and then by presenting results for cylindrical 
panels made of aluminum and laminated-composite 
materials. In particular, the results address the behavior 
of simply supported cylindrical panels that arc subject- 
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ed to combinations of axial compression, transverse 
tension or compression, and shear loads. Results for 
flat plates and curved panels are included. 

Finite-Element Model and Analysis 

The panels considered in this study were analyzed 
with the STAGS finite-element code. 9 The geometry 
and loading conditions and a typical finite-element mod- 
el of a curved panel are shown in Fig. 1. Points on the 
panel mid-surface are located by an x-y-z coordinate 
frame whose origin is at the lower left corner of the pan- 
el. The panel length, width, and radius are defined as L, 
W and R respectively. The panel length and width are L 
= W = 1 4.0 in., for all panels. Flat plates ( R = <*) and a 
shallow curved panel with a radius of R = 60 in. are con- 
sidered in the present study. Simply supported boundary 
conditions are applied to the edges of the panels. The 
panels are loaded by uniform applied traction loads. 
Aluminum and laminated composite panels are consid- 
ered in the investigation. The aluminum panels have a 
nominal thickness of t = 0.12 in. The nominal material 
properties are as follows: Young’s modulus E = 10.0 
Msi, and Poisson’s ratio v = 0.33. The composite pan- 
els are 24-ply laminates. The lamina plies have a nomi- 
nal thickness equal to 0.005 in., which corresponds to a 
nominal panel thickness of t = 0. 12 in. The nominal uni- 
directional properties of a typical 0.005-in-thick lamina 
ply are as follows: longitudinal modulus E j = 19.5 Msi, 
transverse modulus E 2 = 1.45 Msi, in-plane shear modu- 
lus Gj 2 = 0.813 Msi, and major Poisson’s ratio = 
0.30. Three laminate stacking sequences are considered 
in the present study and include a quasi-isotropic [±45/0/ 
90 ] 3s laminate, an axially stiff [±45/0 2 ]3 S laminate, a nd a 
transversely stiff [±45/902]3 S laminate. The finite-ele- 
ment models consist of four-node STAGS 4 1 0 quadrilat- 
eral shell elements. The STAGS 410 element is a flat 
facet-type element that is based on the Kirchoff-Love 
shell theory and the nonlinear Lagrangian-strain tensor. 
The element nodes include three rotational degrees of 
freedom and three translational degrees of freedom. 
Large rotations are accounted for by the use of a co-rota- 
tional algorithm. 

The STAGS code uses both the modified and full 
Newton methods for its nonlinear solution algorithms. 
The equivalence transformation (ET) bifurcation proces- 
sor is used to continue beyond bifurcation points into the 
postbuckling range of the solution. The ET processor al- 
lows the user to select a solution branch in the direction 
of growth of a bifurcation mode that has been computed 
for the nonlinear stress state in the neighborhood of the 
bifurcation point. This bifurcation mode is held constant 
with a magnitude specified by the user, and the load fac- 
tor in the solution is adjusted until the solution residual 


vanishes. 

The prebuckling responses were determined using 
the geometrically nonlinear quasi-static analysis capabil- 
ity in STAGS. The ET processor was used to obtain so- 
lutions beyond the bifurcation point and to branch onto 
the postbuckling equilibrium path of the panels. Once 
the initial postbuckling path was established, the analysis 
was continued into the postbuckling range by using the 
quasi-static analysis capability in STAGS. Linear bifur- 
cation analysis were also conducted using the STAGS 
code. A typical finite-element model contained approx- 
imately 10,000 degrees of freedom. 

Results and Discussion 

The buckling and nonlinear response of panels sub- 
jected to combined loads are presented. The results are 
presented to illustrate the nonlinear behavior of curved 
panels subjected to combined loads. First, results illus- 
trating the response of flat aluminum plates subjected to 
combined loads are presented. Then results illustrating 
the effects of panel curvature on the nonlinear response 
are presented. Finally, results illustrating the effects of 
laminate orthotropy and anisotropy on the buckling loads 
of the panels subjected to combined loads are presented. 

Aluminum Flat Plates 

Nonlinear response curves for geometrically per- 
fect, flat plates subjected to axial compression and trans- 
verse tension or compression loads are presented in Fig. 
2. In particular, five curves are presented in Fig. 2a that 
show the nondimensional load as a function of the aver- 
age end shortening (referred to herein as load-end-short- 
ening curves for convenience) that correspond to the load 
ratio N y /N x = -1, -0.5, 0, 0.5, and I. Similarly, five cor- 
responding curves are presented in Fig. 2b that show the 
nondimensional load as a function of the out-of-plane or 
normal displacement at the center of the plate. These 
curves are referred to herein as load-normal-displace- 
ment curves for convenience. The ordinate is nondimen- 
sionalized by the value of the unidirectional compression 
stress resultant (i.e. N y /N x = 0) in the plate at the bifurca- 
tion point; that is, (N x ) 0 bi f = 325.5 Ibf/in, The end-short- 
ening and panel normal displacements are 
nondimensionalized by the panel length L =14 in. and 
thickness t = 0.12 in., respectively. Negative and posi- 
tive values of the load ratio correspond to transverse ten- 
sion and compression loads, respectively. The filled 
circle on each curve represents the bifurcation point. 

The results in Fig. 2a indicate that the bifurcation 
point of the plate considered increases with increasing 
transverse tension load and decreases with increasing 
transverse compression load, as expected. The values 
of the nondimensional buckling loads range from 0.5 to 
2.08. In addition, the results show a decrease in preb- 
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uckling axial stiffness with increasing transverse ten- 
sion load and the opposite trend is shown for an 
increase in transverse compression load. The results in 
Fig. 2b indicate that the magnitude of the postbuckling 
normal displacement at the center of the plate increases 
with an increase in the transverse compression load and 
decreases with and increase in the transverse tension 
load. The displacement response is characterized by a 
monotonically increasing center displacement with in- 
creasing applied load after buckling. In addition, the 
shapes of the curves indicate that the stiffness of the 
panel increases with decreasing N y (tension). More- 
over, the results in Figs. 2a and 2b indicate that the 
buckling loads obtained from a nonlinear analysis (re- 
ferred to herein as nonlinear-buckling loads) are identi- 
cal to the linear bifurcation buckling loads for all load 
ratios, and these results are summarized in Table 1. 

Nonlinear response curves for geometrically per- 
fect flat plates subjected to axial compression and in- 
plane shear loads are presented in Fig. 3. In particular, 
three load-end-shortening response curves and three 
load-normal-displacements response curves are shown 
for values of the load ratio N xy /N x = ±1, ±0.5, and 0 in 
Figs. 3a and 3b, respectively. The ordinate is nondimen- 
sionalized by the value of the uniaxial compression stress 
resultant (i.e. N xy /N x = 0) in the plate at the bifurcation 
“" point; that is, (N x )°bif* The end-shortening and panel 
normal displacements are nondimensionalized by the 
panel length L and thickness r, respectively. Positive and 
negative values of the load ratio correspond to positive 
and negative in-plane shear loads, respectively. The 
filled circle on each curve represents the bifurcation 
point. 

The results in Fig. 3a indicate that the bifurcation 
point of the plate considered decreases with increasing 
in-plane shear load, as expected. The normalized buck- 
ling loads are 0.86, 0.96 and 1 .0 for load ratios N xy /N x 
= ±1, ±0.5, and 0, respectively. The decrease in buck- 
ling load is attributed to the diagonal compression gen- 
erated by the shear load, however, the shear load has a 
smaller effect on the response than an applied trans- 
verse compression load N y because it produces a small- 
er-magnitude compression load than the transverse 
compression load. The results also show that the preb- 
uckling axial stiffness is the same for all values of in- 
plane shear loads. The results in Fig. 3b indicate that 
the magnitude of the normal displacement at the center 
of the plate decreases with an increase in the in-plane 
shear load. The displacement response is characterized 
by a monotonically increasing center displacement with 
increasing applied load after buckling occurs. Howev- 
er, unlike the results for the transverse compression- 
loaded plates, the postbuckling stiffness does not vary 


substantially with the amount of shear load applied to 
the panel. In addition, the results indicate that the non- 
linear buckling loads for all load ratios agree well with 
the linear bifurcation buckling loads and are summa- 
rized in Table 2. 

Typical bifurcation mode shapes for square alumi- 
num flate plates subjected to combined axial compres- 
sion and transverse tension or compression loads or 
axial compression and in-plane shear loads are shown 
in Fig. 4. The solid and dashed contour lines in the Fig- 
ure indicate positive and negative normal displace- 
ments, respectively and the density of the lines 
indicates the severity of the displacement gradient The 
results indicate that, for panels subjected to uniaxial and 
biaxial compression loads, that is, values of the load ra- 
tio N y /N x = -1, -0.5, and 0, the mode shape is charac- 
terized by a single half-wave response, as shown in Fig. 
4a. In contrast, the results indicate that the mode shape 
is characterized by a single half-wave across the width 
of the plate and two half-waves along the panel length, 
for values of the load ratio N y /N x = 1 , and 0.5, as shown 
in Fig. 4b. A typical bifurcation mode shape for a panel 
subjected to combined axial compression and in-plane 
shear loads (N xy /N x = 1) is shown in Fig. 4c, and is 
characterized by a skewed single half-wave mode shape 
which exhibits a central point of inversion symmetry 
(polar symmetry). The buckling displacements predict- 
ed from a nonlinear analysis of the aluminum flat plates 
are identical to the predicted bifurcation-buckling mode 
shapes, for all loading cases considered. 

Aluminum Curved Panels 

Nonlinear response curves are shown in Fig. 5 for 
the curved panels with R/t = 500 (R = 60 in.) subjected 
to combined axial compression and transverse tension or 
compression. In this figure, six load-end-shortening re- 
sponse curves and six load-normal-displacement re- 
sponse curves are shown for values of the load ratio Ny/ 
N x = -1, -0.5, -0.1, 0, 0.5 and 1. The ordinate is nondi- 
mensionalized by the value of the uniaxial compression 
stress resultant N y /N x = 0 (i.e. N y /N x = 0) in the plate at 
the bifurcation point; that is, (N x )° bif = 723.4 lb/in. The 
end-shortening and panel normal displacements are non- 
dimensionalized by the panel length L =14 in. and thick- 
ness t = 0.12 in., respectively. Negative and positive 
values of the load ratio correspond to transverse tension 
and compression loads, respectively. The open circle 
and square symbols on each curve indicates the bifurca- 
tion point and limit point, respectively. 

The results in Fig. 5a indicate that only the curved 
panel loaded by uniaxial compression, that is N y /N x = 
0, exhibits a bifurcation point. The nonlinear-buckling- 
load value obtained from the nonlinear analysis (N x = 
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720.3 Ib/in.) is within 1% of the load value obtained 
from the linear-bifurcation analysis. Moreover, the 
curved panel with N y /N x = -0. 1 exhibits a limit point at 
a load level of 383.5 lb/in. The limit point for this panel 
is significantly less than the corresponding bifurcation- 
buckling load equal to 727.2 Ib./in. The remainder of 
the panels exhibit a benign, monotonically increasing 
nonlinear response. The normalized nonlinear normal 
displacement response curves for the panels are shown 
in Fig. 5b. The results indicate that the panel subjected 
to uniaxial compression (N y /N x = 0) does not exhibit 
any prebuckling normal deformations. After buckling, 
however, the panel exhibits significant inward radial 
deformations as indicated by the negative values of the 
normalized center displacement response. In contrast, 
the panels subjected to combined axial compression 
and transverse tension or compression do not exhibit a 
bifurcation type response, rather, they exhibit monoton- 
ically increasing load versus normal displacement re- 
sponses from the onset of loading. The panels 
subjected to biaxial compression exhibit positive out- 
ward normal deformations, and the panel subjected to 
combined axial compression and transverse tension ex- 
hibit negative inward normal deformations that have 
the characteristics of initial geometric imperfections. 
The nonlinear, monotonically increasing load-end- 
shortening behavior of the panels, shown in Fig 5a, is 
attributed to the normal deformation response exhibited 
by the panels. More specifically, as the panels deform 
during loading, the geometric stiffness of the panel 
changes which causes the effective axial stiffness of the 
panel to change in a gradual nonlinear manner as indi- 
cated by the nonlinear response shown in Fig. 5a. This 
gradual change in axial stiffness is in contrast to the rap- 
id change in axial stiffness of a panel that exhibits a bi- 
furcation-buckling type response. Overall, these high- 
fidelity results verify the approximate results presented 
in Ref. 8 and indicate generally that linear-bifurcation 
analyses may be misleading for curved panels subjected 
to combined axial compression and transverse tension 
or compression loads. In addition, these results suggest 
that curved panels subjected to this type of combined 
loading condition may not have a high degree of imper- 
fection sensitivity. A summary of the predicted bifur- 
cation-buckling loads and the nonlinear-buckling loads 
for the curved aluminum panels are present in Table 1 . 

Nonlinear response curves for geometrically per- 
fect curved panels subjected to axial compression and in- 
plane shear loads are presented in Fig. 6. Three load- 
end-shortening curves and three load-normal-displace- 
ments curves are shown for values of the load ratio N xy / 
N x = ±1, ±0.5, and 0 in Figs. 6a and 6b, respectively. 
The ordinate is nondimensionalized by the value of the 


uniaxial compression stress resultant (i.e. N xy /N x = 0) in 
the plate at the bifurcation point; that is, (N x ) bif = 723.4 
lbf./in. The end-shortening and panel normal displace- 
ments are nondimensionalized by the panel length L and 
thickness t , respectively. Positive and negative values of 
the load ratio correspond to positive and negative in- 
plane shear loads, respectively. The open circle on each 
curve represents the bifurcation point. 


The results in Fig. 6a indicate that the bifurcation 
point of the compression-loaded panel considered de- 
creases with increasing in-plane shear load, as expect- 
ed. The normalized buckling loads are 0.79, 0.93 and 
1.0 for values of the load ratio N xy /N x = ±1, ±0.5, and 
0, respectively. In particular, the curved panel with a 
load ratio value of N xy /N x = ±1 exhibits a normalized 
buckling load of 0.79 as compared to the corresponding 
flat plate which exhibits a normalized buckling load 
equal to 0.86. This result suggests that panel curvature 
can have a significant effect on the response of the pan- 
el to combined loads. The results also show that the 
prebuckling axial stiffness is the same for all values of 
in-plane shear loads. The bifurcation point is followed 
by an unstable buckling response. In addition, the re- 
sults indicate that the curved panel subjected to com- 
bined axial compression and in-plane shear loads 
exhibit bifurcation-buckling type responses for all load 
ratios considered. In contrast, the corresponding panels 
subjected to combined axial compression and trans- 
verse tension or compression loads exhibit monotoni- 
cally increasing, stable responses. The nonlinear 
normal displacement response curves for the panels are 
shown in Fig. 6b. The results indicate that during the 
prebuckling response the panels do not exhibit any 
normal deformations, however, the post buckling dis- 
placement response of the panels is characterized by 
large magnitude inward normal displacements of the 
central region of the panels. 


A comparison of bifurcation mode shapes and 
buckling deformation patterns predicted from a nonlin- 
ear analysis for square aluminum curved panels sub- 
jected to combined axial compression and transverse 
tension or transverse compression loads or axial com- 
pression and in-plane shear loads are shown in Fig. 7. 
The solid and dashed contour lines in the figure indicate 
positive and negative normal displacements, respec- 
tively, and the density of the lines indicates the severity 
of the displacement gradient. The results indicate that 
the predicted bifurcation mode shape for a panel sub- 
jected to uniaxial compression exhibits a single half- 
wave response and is characterized by a deformation 
pattern ellipse-shaped contours, as shown in Fig. 7a. 
The results also indicate that the buckling deformation 
pattern predicted from a nonlinear analysis of the panel 
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agrees well with the bifurcation mode shape. The re- 
sults indicate that, for a panel subjected to biaxial com- 
pression with a load ratio N y /N x = 1, the panel 

deformation response is characterized by a single half- 
wave response, as shown in Fig. 7b. However, these re- 
sults indicate that there is a slight difference in the de- 
formation pattern determined by a linear bifurcation 
analysis compared to the nonlinear deformation pattern. 
More specifically, the bifurcation buckling mode shape 
is characterized by a deformation pattern with ellipse- 
shaped contours and is in contrast to a deformation pat- 
tern predicted with circular-shaped contours from a 
nonlinear analysis. The results indicate that, for a panel 
subjected to axial compression and transverse tension 
with a load ratio N y /N x = -1, the bifurcation buckling 
mode shape is characterized by two half-waves across 
the width of the plate and one half-wave along the panel 
length as shown in Fig. 7c. However, the predicted de- 
formation pattern from a nonlinear analysis is charac- 
terized by a single half-wave response. A typical 
bifurcation mode shape for a panel subjected to com- 
bined axial compression and in-plane shear loads with 
a load ratio N y /N x = 1 is shown in Fig. 7d, and is char- 
acterized by a skewed single-half wave mode shape 
which exhibits a central point of inversion symmetry. 
The corresponding predicted buckling displacements of 
the panel from a nonlinear analysis agree well with the 
predicted bifurcation buckling mode shape. The results 
in Fig. 7 indicate that the predicted bifurcation buckling 
mode shapes agree with buckling deformation patterns 
predicted from a nonlinear analysis when the nonlinear 
response exhibits a bifurcation response, e.g., panels 
subjected to uniaxial compression loads or combined 
axial compression and in-plane shear loads, and do not 
exhibit a prebuckling normal deformations. However, 
for curved panels subjected to combined loads that do 
not exhibit a bifurcation-buckling type response, the 
predicted bifurcation mode shapes tend to differ from 
the deformation patterns predicted from a nonlinear 
analysis and may be misleading. 

Effects of Laminate Orthotropv and Anisotropy 

Results for three laminates including a quasi-isotro- 
pic [±45/0/90]3 S laminate, an axially stiff [±45/0 2 ]3 S lam- 
inate, and a transversely stiff [±45/90 2 ]3 s laminate were 
obtained to illustrate the effects of laminate orthotropy 
and anisotropy on the nonlinear and buckling response of 
the panels. Similar results to those presented for the alu- 
minum panels were obtained and indicate that, in gener- 
al, the nonlinear and buckling response characteristics 
for the composite laminated panels are similar to those 
exhibited by the aluminum panels. Results illustrating 
the effects of orthotropy on the buckling loads of the flat 


plates and curved panels subjected to combined loads are 
presented in this section. Results for flat plates and 
curved panels subjected to combined axial compression 
and transverse tension or compression loads are present- 
ed in Figs. 8a and b, respectively. Bifurcation-buckling 
loads and buckling loads predicted from a nonlinear 
analysis are denoted by the square and triangle symbols, 
respectively. The bifurcation-buckling values of the load 
N x are denoted by (N x ) bif in the figure, and (N x )° bjf cor- 
responds to the buckling load value for a given panel 
with N y /N x = 0. Similarly, the buckling load predicted 
from a nonlinear analysis are denoted by (N x ) cr in the 
figure and (N x )° cr corresponds to the buckling value for 
a given panel with N y /N x = 0. 

The results in Fig. 8 show that panel curvature and 
laminate orthotropy can have a significant effect on 
buckling response of the panels subjected to combined 
axial compression and transverse tension loads (load ra- 
tio N y /N x = -1, and -0.5). More specifically, the results 
in Fig. 8a indicate that axially stiff composite flat plates 
have higher values of normalized buckling loads that 
range from 1.0 to 2.46 followed by the results for the 
aluminum and quasi-isotropic composite panels, which 
exhibit nearly identical nondimensional buckling loads, 
with normalized buckling load values that range from 
1 .0 to 2.06 for the loading conditions considered, and fi- 
nally the results for the transversely stiff composite flat 
plates with normalized buckling load values that range 
from 1.0 to 1.52. In contrast, the results indicate that 
the flat plates exhibit the same nondimensional buck- 
ling loads when subjected to uniaxial and biaxial com- 
pression loads (load ratio N y /N x = 0, 0.5, and 1). In 
addition, the results indicate that the nonlinear-bucking 
loads agree well with the predicted bifurcation-buck- 
ling loads for these flat plates and the values are gener- 
ally with in 1%. The results presented in Fig 8b indicate 
that the curved panels with R/t = 500 exhibit a similar 
trend; that is, the nondimensional buckling loads are 
generally higher for the axially stiff composite panels 
followed by the results for the aluminum and quasi-iso- 
tropic composite panels and finally the results for the 
transversely stiff composite panels. However, the de- 
gree of variation in the values of the normalized bifur- 
cation-buckling loads is much less than that exhibited 
by the fiat plates. These results suggests that a nonlin- 
ear coupling exists between the panel curvature and 
laminate orthotropy for the panels considered. In addi- 
tion the results indicate that the curved panels subjected 
to the combined axial compression and transverse ten- 
sion or compression loads do not exhibit a bifurcation 
point, thus these bifurcation-buckling results may be 
misleading since bifurcation buckling may have no 
meaning in the context of curved panels subjected to 
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these types of combined loads. 

Results are presented for flat plates and curved 
panels subjected to combined axial compression and in- 
plane shear loads in Figs. 9a and b, respectively. The 
results indicate that panel curvature and laminate 
orthotropy and anisotropy can have an effect on the 
nondimensional buckling load values for the loading 
conditions considered. More specifically, the results in 
Fig. 9a suggest that, for the flat plates, panel orthotropy 
has little effect on the nondimensional buckling load 
level as is evidenced by that fact that all the results for 
the composite flat plates are nearly identical for all val- 
ues of the load ratio N xy /N x . However, the results also 
suggest that the response is affected by laminate anisot- 
ropy. More specifically, the results show that the buck- 
ling load value is dependent on the sign of the load 
ratio, that is the buckling load value is dependent on 
whether the in-plane shear load is positive or negative. 
For example, The composite plates with a load ratio 
N xy /N x = -1, corresponding to a negative in-plane shear 
load, exhibit a nondimensional buckling load of 0.835, 
however, the composite plates with a load ratio N xy / 
N x = +1, which corresponds to a positive in-plane shear 
load, exhibit a nondimensional buckling load of 0.875. 
This result is in contrast to the results for the aluminum 
plate in which the buckling loads are independent of the 
sign of the load ratio. The results for the curved panels 
subjected to combined axial compression and in-plane 
shear loads are presented in Fig. 9b. The result indicate 
a similar trend to that exhibited by the corresponding 
flat plates. However, the results indicate that, in gener- 
al, the curved panels exhibit lower normalized buckling 
loads when subjected to combined loads than do the 
corresponding flat plates. In addition, the results indi- 
cate that laminate anisotropy affects the buckling loads 
of the curved panels in a similar manner to that exhibit- 
ed by the flat plates as evidenced by the dependence of 
the buckling load on the sign of the in-plane shear load. 
Moreover, the results indicate that laminate orthotropy 
of the panels considered also has a slight effect on the 
buckling response. More specifically, the results indi- 
cate that the transversely stiff panels have higher values 
of normalized buckling loads which range in value from 

0.8 to 1.0 followed by the axially stiff, quasi-isotropic, 
and aluminum panels with values of normalized buck- 
ling loads that range from 0.77 to 1.0. 

Concluding Remarks 

High-fidelity analysis tools and advanced compu- 
tational capability now make it possible to develop reli- 
able, high-fidelity design criteria for the buckling 
behavior of shell structures subjected to combined me- 
chanical and thermal loads. As part of updating and de- 
veloping such design criteria, in-depth examination of 


the nonlinear response of shells subjected to combined 
loads is needed. The present paper addressed this basic 
need by presenting results and discussion for flat plates 
and curved panels subjected to a variety of combined 
loading conditions. In particular, results for loading 
conditions that involve combinations of axial compres- 
sion, transverse tension or compression, and shear loads 
have been presented. Results that illustrate the effects 
of laminate orthotropy and anisotropy have also been 
presented. The results indicated that panel curvature 
can have a significant effect on the response of panels 
subjected to combined loads. In particular, curved pan- 
els subjected to combined axial compression and trans- 
verse tension or compression loads do not exhibit a 
bifurcation-buckling type response. As a result, a linear 
bifurcation analysis could misrepresent the physical be- 
havior of the panel. The results also indicated that lam- 
inate orthotropy and anisotropy can have a substantial 
effect the buckling loads of the flat plates and curved 
panels subjected to combined loads. These results 
should be of interest to engineers designing practical 
shell structures. 
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Table 1: Buckling loads for aluminum panels subjected to combined axial compression and transverse 

tension or compression loads. 


Loading 

Ratio, 

N/N x 

R/t = co 

(Flat plate) 

R/t = 

500 

Bifurcation- 
buckling 
l03d (N x ) bif , 
lbf/in. 

Nonlinear- 
buckling 
load (N x ) cp 
lbf/in. 

Bifurcation- 
buckling 
load (N x ) bify 
lbf/in. 

Nonlinear- 
buckling 
load (N x ) cp 
lbf/in. 

-1.0 

677.7 

677.7 

1215 

a 

-0.5 

581.1 

581.1 

1051 

a 

-0.1 

- 

- 

727.2 

383. 5 b 

0.0 

325.6 

325.6 

723.4 

720.3 

0.5 

217.0 

217.0 

496.5 

a 

1.0 

162.8 

162.8 

376.3 

a 


* Monotonically increasing response which does not exhibit a bucking point 
b Limit point 


Table 2: Buckling loads for aluminum panels subjected to combined axial compression and in-plane shear 

loads. 


Loading 

Ratio, 

N xy /N x 

R/t = oo 

(Flat plate) 

R/t = 

500 

Bifurcation- 
buckling 
load (N x ) bif , 
lbf/in. 

Nonlinear- 
buckling 
load (N x ) cp 
lbf/in. 

Bifurcation- 
buckling 
load (N x ) bif > 
lbf/in. 

Nonlinear- 
buckling 
load (N X ) CP 
lbf/in. 

-1.0 

281.1 

281.1 

570.9 

568.5 

-0.5 

312.0 

312.0 

674.8 

666.1 

0.0 

325.6 

325.6 

723.4 

720.3 

0.5 

312.0 

312.0 

674.8 

666.1 

1.0 

281.1 

281.1 

570.9 

568.5 
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Fig. 1 Typical panel geometry and loading conditions and finite-element model. 



A/L 5 ft 

a) Load-end-shortening response b) Load-normal-displacement response 


Fig. 2 Nonlinear response of square aluminum flate plates subjected to combined axial compression and transverse ten- 
sion or compression loads (L = 14.0 in., r = 0*12 in., (N x )° bif = 325.5 lbf/in.). 
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Fig. 3 Nonlinear response of square aluminum flate plates subjected to combined axial compression and in-plane shear 
loads (L = 14.0 in., t = 0.12 in., (N x )° bif = 325.5 Ibf/in.). 



b) Ny/N x = -1, -0.5 




Fig. 4 Topical bifurcation mode shapes for square aluminum flat plates subjected to combined axial compression and 
transverse tension or in-plane shear loads. 
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a/l m 

a) Load-end-shortening response b) Load-normal-displacement response 


Fig. 5 Nonlinear response of square aluminum curved panel subjected to combined axial compression and transverse ten- 
sion or compression loads (L = 14.0 in., R = 60 in., t = 0.12 in., (N x )°bif = 723.4 Ibf/in.). 



A/L 5 It 

a) Load-end-shortening response b) Load-normal-displacement response 

Fig. 6 Nonlinear response of square aluminum curved panel subjected to combined axial compression and in-plane shear 
loads (, L = 14.0 in., R = 60 in., t = 0.12 in., (N x )° bir = 723.4 lbtfin.). 
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■ 


■ 


bifurcation mode shape buckling deformation pattern 
a) Uniaxial compression, Ny/N x = 0 

■ 


B 




bifurcation mode shape buckling deformation pattern 
b) Biaxial compression, N y /N x = 1.0 





bifurcation mode shape buckling deformation pattern bifurcation mode shape buckling deformation pattern 

c) Axial compression and transverse tension, Ny/N x = -1 d) Axial compression and in-plane shear, N xy /N x = 1 

Fig. 7 Comparison between bifurcation buckling mode shapes and buckling deformation patterns predicted from a non- 
linear analysis for curved aluminum panels subjected to combined axial compression and transverse tension or compres- 
sion loads or combined axial compression and in-plane shear loads. 


□ Bifurcation buckling load 
a Nonlinear buckling load 




Fig. 8 Effects of orthotropy and curvature on the nondimensional buckling loads for square panels subjected to combined 
axial compression and transverse tension or compression loads. 
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Fig. 9 Effects of orthotropy and curvature on the nondimensional buckling loads for square panels subjected to combined 
axial compression and in-plane shear loads. 
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